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Abstract 

Gribov's scenario of supercritical charges in QCD is investigated. We perform 
a numerical study of the corresponding equation for the Green function of light 
quarks. This is done in an approximation which neglects all pion contributions. 
Different types of solutions in the Euclidean region are discussed and the mass 
function of the quark is calculated. The solutions of the equation are shown to 
have a qualitatively different behaviour if the strong coupling constant exceeds 
a critical value ac — 0.43 in the infrared region. Chiral symmetry breaking is 
found to occur at supercritical coupling. The analytic structure of the solutions is 
investigated. Earlier results obtained by Gribov are confirmed and extended. 
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1 Introduction 



The breaking of chiral symmetry and the confinement of quarks and gluons are two of 
the most important properties of QCD. The details of the mechanism leading to confine- 
ment are still largely unknown, and the understanding of non-perturbative dynamics 
in QCD in general is still rather poor. A new picture of the confinement mechanism 
and of chiral symmetry breaking was developed by V.N. Gribov It is based on 

the phenomenon of supercritical charges which can occur in QCD due to the existence 
of very light quarks. Its consequence is a dramatic change in the vacuum structure of 
the light quarks compared to the usual perturbative picture at small coupling. 

The phenomenon of supercritical charges is well-known in QED (for an extensive 
review see ||5|). The energy of the bound-state levels in the field of an isolated heavy 
nucleus decreases if the charge Z of the nucleus is increased. When the charge exceeds 
a critical valueQ of Z^r = 137, the lowest bound-state level dives into the Dirac sea, i. e. 
sinks below —rrie- As a consequence an electron from the (filled) continuum undergoes 
a transition into this level, and a positron is emitted. The electron is said to 'fall onto 
the center'. In this situation the simple quantum mechanical picture breaks down, and 
the emerging bound state is in fact a collective state with a high probability to find 
an electron very close to the nucleus. This mechanism is called supercritical binding. 
The condition for its occurrence is that the Compton wavelength 1/m of the electron 
is much larger than the radius of the heavy charge. 

Gribov's confinement scenario is based on the idea that a similar phenomenon occurs 
in QCD due to the existence of very light (almost massless) quarks. The crucial point 
is that in this scenario already the color charge of a single quark is supercritical. Since 
this applies also to the light quarks themselves the situation is more involved than in 
QED. We will give only a condensed description of the resulting scenario here. More 
detailed accounts have been given in [§, ^ 0. In order to get an understanding of the 
underlying physical picture we again use the quantum mechanical description, having 
in mind that the quantitative analysis should be based on the full underlying quantum 
field theory. 

The confinement of heavy quarks in Gribov's scenario is very similar to the super- 
critical binding in QED. Due to its supercritical charge the heavy quark captures a 
light antiquark from the vacuum, thus decaying into a supercritical heavy-light bound 
state. At the same time a light quark is created. This light quark decays again, as we 
will discuss now. 

In order to understand the confinement of light quarks, we first consider a bound 
state of a light quark and a light antiquark. If the coupling constant is small this is just 
a normal bound state like positronium in which the quarks have positive kinetic energy. 
If we now increase the coupling the binding energy will also increase. The total energy 
of this state will thus decrease. But if the coupling is further increased — and becomes 
supercritical — a situation is possible in which the total energy of the bound state 
becomes negative. In order to have a stable vacuum, however, the existence of negative 
energy particles has to be avoided. Consequently, the corresponding quark states in 
the supercritical bound states have to be filled in the vacuum. Therefore there are 
filled quark states with positive kinetic energy in addition to the usual filled states in 

'^This number holds for a point-like nucleus, for an extended charge it is around Zc- ~ 165. 
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the Dirac sea (see^ Fig. |T|). The scale fip separating filled and empty states of positive 




Figure 1: Vacuum of light quarks in Gribov's scenario, and q denote the (kinetic) 
energy and three-momentum, respectively. 

kinetic energy resembles the Fermi surface in solid state physics. The existence of the 
additional states in the vacuum of light quarks implies also the existence of additional 
excitations of this vacuum which have quite unusual properties. The qq pair of such an 
excitation forms a supercritical bound state in which the quark and antiquark both have 
negative kinetic energy. They are interacting repulsively, and the supercritical bound 
state has positive total energy. The 'binding force' leading to this unusual meson^is the 
Pauli exclusion principle. The quark and antiquark are bound in this meson because 
all other energetically possible states in the vacuum are already filled. 

Having discussed the emergence of the novel meson states, we can now understand 
the confinement of light quarks. According to Gribov it is caused by the continuous 
decay of the light quark. Any quark (with positive or negative kinetic energy, or 
(?(_)) decays into a supercritical bound state M and a quark of negative kinetic 

^In this figure only the energy of the quark is shown. One has to keep in mind that some of the 
states shown here exist only within supercritical bound states with an antiquark as described above. 
^Some possible properties of these novel mesons have been discussed in M H. 
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energy, 

q{±) + (?(_) . (1) 

In this sense the quark exists only as a resonance and cannot be observed as a free 
particle. 

Since the existence of the novel meson states is due to the Pauli principle it is 
immediately clear that the above confinement mechanism works only for quarks but 
not for gluons. But the confinement of gluons could possibly be a 'second order effect' 
in Gribov's scenario, namely due to their coupling to light quarks which subsequently 
decay as described above. 

It is obviously desirable to find a quantitative description for this interesting physical 
picture of supercritical color charges. The confinement of quarks and gluons should be 
encoded in the singularities of the respective Green functions. Therefore the Green 
function of the quark is a suitable object to study in this context. In ||2|, |^ Gribov 
derived an equation for the retarded Green function of light quarks. It takes into 
account especially the dynamics of the infrared region but also reproduces asymptotic 
freedom at large momenta. Chiral symmetry breaking has been found to occur when the 
strong coupling constant exceeds a critical value, leading to the emergence of Goldstone 
boson (pions). It has been argued in ^ that the nature of these Goldstone bosons is 
such that they should in fact be regarded as elementary objects. Corrections to the 
Green function caused by these Goldstone bosons are expected Q to lead a Green 
function of light quark which exhibits confinement, whereas the equation without these 
corrections is not expected to imply confinement [P, |^. Unfortunately, the paper Q 
remained unfinished, and a full study of the analytic properties of the Green function 
and their consequences still remains to be done. Especially, it will be important to 
see how and to what extent the — though somewhat simplified — physical picture 
described above can be derived from the analytic properties of the resulting Green 
function. 

Gribov's equations for the Green function of light quarks (with or without pion 
corrections) are non-linear differential equations. So far the studies of these equations 
[^, 1^ have been performed only by means of asymptotic expansions. It is the purpose 
of the present paper to perform a complete numerical study of the equation without pion 
corrections. This allows us to study the breaking of chiral symmetry also quantitatively 
and in more detail. We also investigate the analytic structure of the Green function 
resulting from Gribov's equation. 

It turns out that in Gribov's equation the critical value of the strong coupling 
constant is surprisingly low, Ofc = 0.43. In the derivation of the equation it is assumed 
that the coupling constant does not become very much larger than this critical value. To 
some extent Gribov's approach can thus be considered as a semi-perturbative approach 
to confinement. This picture might also explain why we observe an essentially smooth 
behaviour of non-perturbative effects in the transition from the parton level to the 
hadron level. Typical multiplicities at the parton level, for example, are in surprising 
correspondence to those observed at the hadron level (for a more detailed discussion of 
this and similar observations see The idea of an infrared finite coupling has also 

been widely discussed in the context of power corrections and the dispersive approach 
to renormalons in QCD |11, see also |13| and references therein. In that approach, it 



appears consistent to define an effective running coupling down to very low momentum 
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scales, in the sense that its first moments have a universal meaning. The values of the 
coupling found in the corresponding experimental analyses are in fact bigger than the 
critical value in Gribov's equation. 

The equation can be derived from the Dyson-Schwinger equation for the Green 
function of the quark in Feynman gauge. The approximations made in the derivation 
are motivated by the underlying physical picture, especially concerning the behaviour 
of the strong coupling constant. This method can therefore also be viewed as an uncon- 
ventional approach to the difficult problem of solving the Dyson-Schwinger equations in 
QCD (for a review see 0]). The approximations usually made in solving the Dyson- 
Schwinger equations are intrinsically difficult to control. Comparisons with results 
obtained in Gribov's approach will therefore be potentially very useful. 

The paper is organized as follows. In section ^ we outline the main steps leading 
to the equation for the Green function of light quarks and describe some of its most 
important properties. In section ^ a suitable parametrization of the Green function 
is given. The asymptotic behaviour of the equation for small and large momenta is 
discussed and the critical value of the strong coupling constant is derived. Section |^ 
deals with the Euclidean region of space-like momenta. In section 4.1 the solutions 



are shown to exhibit asymptotic freedom at large space-like momenta. Section 4.2 



provides models for the running coupling at small (space-like) momenta which are 
needed for the numerical analysis of the equation. The possible types of solutions in 



the Euclidean region are classified in section 4.3. The characteristic change in the 



solutions at supercritical coupling is discussed. Section ^ deals with the behaviour of 
the dynamical mass function of the quark in the Euclidean region. Phase transitions 
are found to occur for supercritical coupling and lead to chiral symmetry breaking. 
We study how this effect depends on the models used for the running coupling in the 
infrared. In section |^ we determine the analytic structure of the solutions in the whole 
momentum plane for the different types of solutions classified in section O. We close 
with a summary and an outlook. 

The results presented in sections 3^ and iA concerning the asymptotic behaviour 
of the equation have partly been obtained already in They have been included in 
some detail in the present paper since they are immediately relevant to our analysis 
and serve to make it self-contained. 



2 The equation for the Green function of Ught quarks 

In this section we will outline the main steps that lead to the equation for the Green 
function of light quarks and and highlight some of its properties which are relevant to 
our discussion. Some of these properties and the full derivation of the equation have 
been discussed in detail in [§, ^. 

The first step is the choice of a gauge. As noted by Gribov, the Feynman gauge 
turns out to be particularly well suited for deriving a simple equation which is especially 
sensitive to the infrared dynamics. In other gauges it would be extremely difficult to 
find a similarly simple equation. The physical results, like for example the occurrence 
of chiral symmetry breaking, will of course be independent of the choice of gauge. In 
Feynman gauge the gluon propagator has the form 

D,Ak) = -^as{k')- (2) 
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The exact behaviour of the strong couphng constant is not known at small momenta. 
In Gribov's derivation of the equation it is assumed that the coupling constant is a 
slowly varying function of the momentum and does not become very large at small 
momenta. Such a behaviour is sketched in Fig. It turns out that the occurrence 




Figure 2: Assumed behaviour of the strong coupling as ((7^) 

of a supercritical behaviour of the Green function does not depend on the details of 
the coupling in the infrared, as long as its value is above the critical value in some 
interval of momenta. As we will see, this critical value is rather low, Uc = 0.43. 
These properties of the running coupling are consistent with the picture arising in the 



dispersive approach ||ll|] to power corrections in QCD (for reviews see |10, |l^]). There 
it appears that the definition of a running coupling constant at very low momenta is 
possible in the sense that its integral moments have a universal meaning. Motivated by 



this, possible models for the coupling have been constructed, see for example ||15|, [L6 |. 
For our numerical study we will choose a rather simple form of the coupling, see section 
L% below. 

One starts from the Dyson-Schwinger equation for the inverse Green function G^^ 
of the quark and considers its perturbative or diagrammatic expansion. To the cor- 
responding sum of diagrams one applies the double differentiation = d^d^, where 

is the derivative with respect to the external momentum of the quark. Firstly, 
this is a way to regularize the divergences in these diagrams, and gives a finite result. 
Secondly, it can be used to collect the most singular contributions to the quark Green 
function from the infrared region. This is based on the observation that the action of 

on the gluon propagator in Feynman gauge gives a delta function, 

1 ....... . 



V\2 



{q-q' 

The integration variables in all diagrams can be arranged in such a way that the external 
momentum of the quark is carried along gluon lines. If a gluon line is now differen- 
tiated twice, the above identity then transforms the integral over the corresponding 
gluon momentum into two zero-momentum gluon insertions. All other contributions 
to the respective integrals — those with derivatives in different lines or of the running 
coupling — are clearly less singular in the infrared region. It is in principle possible 
to treat these terms systematically as corrections. But the resulting equation will then 
be a complicated integro-differential equation, as briefly indicated in |§. In first ap- 
proximation those contributions are neglected, and one is left with a series of diagrams 
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with two gluon insertions that carry zero momentum. This sum can be shown to be 
the diagrammatic expansion of a full inverse quark Green function with two full quark- 
gluon vertices r^(g, /c = 0) inserted. Using Ward identities the latter can be replaced 
by derivatives d^G~^ of the quark Green function. Having eliminated the vertex func- 
tions, one ends up with a second order differential equation for the Green function of 
a light quark, 

d^G-' =g{d''G-')G{d^G-'), (4) 

where 

g = CF • (5) 

TT 

This is Gribov's equation which will be the subject of our study. 

A comment is in order concerning the choice of scale of the running coupling in 
@. As described so far, the derivation of the equation has concentrated on the most 
important contributions from the infrared region. But it is of course desirable to find 
an equation which describes the Green function correctly also in the ultraviolet region. 
The use of the relation implies that the coupling has to be evaluated at zero momen- 
tum. But it can be shown that by replacing 0^(0) by as{q) one arrives at an equation 
that also reproduces the correct behaviour at large momenta. Given the assumptions 
about the running coupling discussed earlier, the correction induced by this replace- 
ment is sub leading as far as its contribution to the infrared region is concerned. In the 
approximation presently considered we can therefore accept equation (^) with (^) as 
an equation that is expected to provide an adequate description of the Green function 
at all momentum scales. 

For simplicity, equation (Q) is written for one-flavour QCD. This is sufficient as 
long as we are mainly interested in the occurrence of confinement and chiral symmetry 
breaking. The generalization to the more realistic case of a doublet of light quarks is 
straightforward and will be important for the study of bound states in Gribov's picture. 

The fact that equation (Q) is a second order differential equation implies that its 
solutions will involve two dimensionful constants of integration. These will be related 
to the quark mass and the quark condensate. 

An obvious and important property of the equation is its invariance under a rescaling 
of the Green function, G cG for any constant c. As a consequence, the equation will 
not involve the full wave function renormalization but only its logarithmic derivative. 
The equation is not scale invariant with respect to the momentum. The breaking of 
scale invariance is due to the presence of the running coupling constant. It is only 
through the running of the coupling that a momentum scale is introduced. 

A further property of the equation is a certain symmetry between the Green function 
G and its inverse G~^. One can easily show that the equation (Q) implies 

d^G = i2-g){d^G)G-\d^G). (6) 

This means that G solves the same equation as G~^, but with g replaced hy 2 — g, the 
symmetry point being g = 1. At g = 2, corresponding to Og = 37r/2, the Green function 
would thus become a free Green function. This symmetry is certainly unphysical, 
and we should trust the equation only for comparatively small values of the coupling, 
roughly speaking below one or two. This is in agreement with the fact that in the 
derivation of the equation the coupling was assumed to be small. 
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3 Parametrization and asymptotic behaviour 
3.1 Parametrization 

Due to invariance under parity and Lorentz transformations the inverse Green function 
has the general form 

G-\q)=a{q^)4-h{q^) (7) 
with two scalar functions a and b. We will in the following use the variable 



q = VQ^^ (8) 

such that the half plane 5fteg > already covers the full plane in q'^, the variable in 
which the Green function is usually discussed. It will be convenient to use instead of 
(0) the following parametrization of the Green function^, 

G~' = -p exp (-14) (9) 



2 g 

with two complex functions p and (j). This corresponds to 

aiq ) = - /3 sinh — (10) 
q 2 

h{q^) = /J cosh ^ (11) 

in the parametrization (^). The dynamical mass function M of the quark is then given 
by the function only, 

M(g2) = ^ = ,coth|, (12) 

«(r) 2 

whereas the function p represents the wave function renormalization. In terms of the 
usual notation we have = p/q. We further introduce 

i = \nq = \n./(fq'^ (13) 

and denote the derivative with respect to this variable as 

f{q) = d^fiq) , (14) 

Since the solutions of Gribov's equation (^) depend only on the logarithmic derivative 
of the wave function renormalization, it is useful to define 



where 



p = l + P-. (15) 
P 

(3=l-g = l-CF—. (16) 
vr 

Then the equation ^ for the Green function translates into a pair of coupled differential 
equations for p and cj), 

p=l_p2_^2g^2^3gij^l^20 (17) 

4> + 2p^-3smli4> = 0, (18) 
which will be the basis of our further analysis. 

''This parametrization deviates from the one used in |H, W- 
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3.2 Asymptotic behaviour 

We now study the asymptotic behaviour of the solutions of Gribov's equation. An 
important outcome of this study will be the determination of the critical coupling at 
which chiral symmetry breaking occurs. 

First we keep the coupling constant fixed. The running of the coupling can then be 
treated under the assumption that the asymptotic behaviour of the solutions depends 
smoothly on the coupling. This assumption will be justified by our numerical analysis 
further below. Since the equations ([l^) and ( [l8| ) depend only on the logarithm of q 
the equations are the same along all straight lines passing through the origin of the 
complex g-plane. The initial conditions, at q = for example, do not exhibit this 
apparent symmetry such that the solutions will be different in different directions in 
the g-plane. The fixed points of the equation, however, turn out to be independent of 
the direction in the g-plane. 

Behaviour for |g| — >^oo 

For large \q\ the pair of equations (0), (|l8|) has stable fixed points at 



4> = {2n + l)iTT (neZ); p= ^1 + 3/32. (19) 



As we will see in section 4.1 the existence of these fixed points implies the asymptotic 
freedom of the corresponding solutions. 

The periodicity of the above fixed points is obvious from the equations, and we will 
now concentrate on the fixed point at (j) = in. Perturbing the solutions around the 
fixed point, 

(f) = in + il; ■ p = pQ-\-p (20) 
and expanding to first order in the perturbations we find 

= 1 + 3/?2 > . (21) 

Further we have 

d^p=-2pop, (22) 
such that p = L>e~^^°^ with some D G C For a stable fixed point we thus have to 



choose the positive root po = \/l + 3/3^ > in (pi|). We note that the function p 
consequently develops a singularity 



(23) 



p ~ exp 

The linearized equation for ■(/; becomes 

+ 2poip + 3ip = 0. (24) 
Thus -0 = C7ie^+« + C2e^-« with Ci, C2 G C. We find 



7± = -PO ± ^J3f3^ - 2 . (25) 

Here 7+ and 7_ can be real (/3^ > 2/3 ) or complex (/J^ < 2/3). As we will see in 
section ^ these two possible cases have quite different physical consequences. 
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In the first case, 0^ > 2/3, the function (p approaches iir monotonically. This case 
is characterized by 

g<g^ = I -^^0.18 (26) 

or 

5>l + y|~ 1.82. (27) 

Here we find the critical value of the coupling constant already mentioned earlier, 

37r 

a, = —5,- 0.43, (28) 

at which the solutions change their behaviour. 

In the second case, < 2/3, the function oscillates while approaching the fixed 
point iir. The corresponding supercritical behaviour of the equation is characterized 
by values of the strong coupling constant a<j in the interval 

ac<a. <4.3. (29) 

The emergence of the upper limit is in agreement with the symmetry of the equation 
relating the Green function to its inverse for g ^ 2 — g, see section ^ The fact that the 
equation exhibits subcritical behaviour at very large values of the coupling is certainly 
unphysical. We cannot expect that the equation describes the Green function correctly 
also at very large coupling. 



Behaviour for |g| — ^0 

For small \q\ there are two possible cases. In the first case, (j) approaches one of the 
fixed points described above, (f) = (2n + l)i7r. To show this one can proceed as in the 
case of large \q\. But now we are considering ^ — > — oo and therefore have to choose 
the negative root in ( [2l] ) in order to find a stable fixed point, po = — ^^1 + 3/3^ < 0. 
Again, the solutions will oscillate while approaching the fixed point if the coupling is 
supercritical. 

The other case possible for \q\ — > is that (j) vanishes at q = 0. Linearizing the 
equation for small (p results in 

4> + 2po4> -3(1) = 0. (30) 

Prom the equation for p we find that for this fixed point p ^ po with Pq = 1. With 
the ansatz cp = Ce"'^ it is required that 7 > for (p to be regular. In order to have a 
solution which at large q approaches ivr with damping we need pq = 1 and therefore 
7 = 1. The easiest way to see this is from eq. ( |3^ ) and the corresponding discussion in 
section ^ below. 



Running coupling 

The running of the coupling can be treated assuming that the asymptotic behaviour of 
the solutions depends smoothly on it. The values of the function p at the fixed points 
discussed above depend on (3. Therefore they are changed accordingly, i. e. have to 
be replaced by f5{q = 0) or I3{q = 00) = 1, respectively. The oscillations occurring 
for |q| ^ cxD stop at the scale at which the coupling becomes smaller than a,, and 
approaches the fixed point monotonically above this scale. 
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4 Solutions in the Euclidean region 



We first investigate the equation in the Euchdean region, i. e. for space-hke momenta. 
Therefore we want to consider purely imaginary values of our variable g, thus q = iq 
with a real-valued and positive q. The derivative with respect to q will be denoted 

= f'iq) ■ (31) 

aq 

For space-like momenta the dynamical mass function M{q^) is required to be real- 
valued. Eq. ( [l^ then implies that cj) is purely imaginary. In addition, a real-valued 
mass function requires that the function p is real-valued for space-like momenta. For 
convenience we define for the use in the present section 

4> = iX. (32) 

where x is real-valued. In this section we thus have to consider only real-valued 
functions x ^'^^ P depending on the real parameter q. 



The equation (18) for x (or (/>, respectively) can be reformulated in such a way that 



it permits a simple interpretation. The function 

e = Y - 3(1 - cosx) (33) 

can be interpreted as the energy of a motion with x being a one-dimensional degree of 
freedom. The equation of motion equivalent to ( p^ ) is 

d^e = -2px^ . (34) 

The behaviour of x can then be interpreted as a motion with damping (given by p) in 
the potential 

y = -3(1 -cosx)- (35) 

This potential has minima at x = (2n + l)7r for all n G Z. Thus for space-like momenta 
the fixed points discussed in section |3.2| appear as the minima of the potential V. 



4.1 Asymptotic freedom 

For large = —q^ the Green function should behave according to perturbative renor- 
malization and exhibit asymptotic freedom. We now show that Gribov's equation 
reproduces exactly this behaviour for solutions that approach one of the fixed points 
discussed above. This is done by considering the leading terms in the limit of large q. 

We first consider the wave function renormalization. The corresponding renormal- 
ization constant is in our parametrization defined as 

p = e^Z-\i)=qZ-\ (36) 



as can be seen when eq. (IC) is evaluated at ~ in. As is usually done we assume Z{^) 
to be a slowly varying function. Using (|l5|) and ([l^) one derives in the limit X ~^ the 
following equation^] for /), 

1 



{d^ + 3)(ag - l)p + 35/9 - g-{d^pY = . (37) 



^In Q this equation (there eq. (4.42)) contains a misprint. 
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Inserting ( ^6| ) and neglecting terms of the order c?|Z ^ and {d^^Z we find 

d,Z~^ + IgZ'^ = 0. 



(38) 



This coincides with the well-known wave function renormalization in Feynman gauge 
as it can be found for example in |17]. Turning to mass renormalization we observe 
that close to one of the fixed points, (p = ix = — o") with small a, the mass function 
is given by (see eq. (p^)) 



a = 2e~«M(e) 



2M 



(39) 



Linearizing eq. (|l^) we find 



We further note that due to (l3q) and (|3q) 



l--g. 



P 

P 



1 + 



(40) 



(41) 



Together with (^) this can be inserted in (|40|). Neglecting the term of order d'jM we 



arrive at 



dfM 



zgM, 



(42) 



which is exactly the mass renormalization at one loop. Its solution is 



M(g2 



mo 



(43) 



where ttiq is the mass at a given scale goi and in one-loop approximation the exponent 
is 7m = 4/6o with = 1\ — |nj. 

Taking into account also the sub-leading solution for a from ( ^0[ ) or, equivalently, 
from (U) we find that a behaves at large q as 



2m„ 

a ~ -\ 

q q-^ 



(44) 



as can be seen from (25) since in this limit j3 ^ 1. Accordingly, the mass function 
behaves as 



Pi'2 



(45) 



We thus find two dimensionful parameters, niq and v^. These can be identified with 
the quark mass and a quark condensate, respectively. It is difficult to disentangle the 
two terms in (|4^) numerically. We will therefore not pursue this interesting issue any 
further in the present paper. 
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4.2 Models for the running coupling 

The behaviour of the strong couphng constant is only very vaguely known in the infrared 
region. In order to perform a numerical study we have to use a model for Og- The model 
is required to be in agreement with the general assumptions used in the derivation of 
the equation (see section |2|). Obviously, any useful model should coincide with the 
perturbative running of the coupling at large momentum scales. In the following we 
will use two different models of this kind. Both are rather simple but should be sufficient 
for studying the physical effects resulting from Gribov's equation. More complicated 



models (see for example [15, |16[]) could of course be implemented in the same way. 

For the perturbative behaviour of the coupling we have to specify Aqcd and Uf. 
Moderate changes in these parameters do not have any significant effect on our results 
for the Green function since they only change the behaviour of at large momenta 
where the coupling is subcritical. To be specific we choose Aqcd = 250 MeV and n/ = 3. 
(The latter choice is not completely in agreement with the fact that the equation is in 
the present paper studied for only one flavour. The choice Uf = 1 would lead to almost 
identical results for the Green function.) 

Our first model is the more realistic one, and is motivated by the dispersive approach 
to power corrections in QCD This approach is based on the assumption that the 
coupling constant can be defined down to very small momenta, and that this coupling 
in the infrared has a universal meaning. Then it is possible to determine its integral 
over the infrared region from measurements of infrared and collinear safe observables 
like certain event shape variables (for a recent review see [|lO|] ). In this way one finds 
for the integral of the coupling 

"°"2Givi) (^s(k)dk^ 0.5. (46) 

This condition can be fulfilled by shifting the argument of the logarithm in the usual 
one-loop formula for the perturbative running of the coupling, 

= 7 • (47) 

(11 - |n/j ln(-(?2/A^cD + a) 

For a = this is exactly the one-loop renormalization of the coupling. Our first model 
for the running coupling is obtained for a = 6. This choice is made to satisfy the 
condition (p6|). We will in the following refer to this model as type A. This running of 
the coupling is shown as curve A in Fig. ^. For comparison we show in that figure as 
curve C also the coupling obtained from the one-loop renormalization. 

The second model is shown as curve B in Fig. |^. It is obtained from the one-loop 
renormalization of the coupling (see eq. (^) with a = 0) by simply cutting it off at some 
given value and assuming it to be constant below the corresponding momentum scale. 
In order to avoid problems in the numerical treatment of the equation we smooth out 
the resulting edge by fitting a polynomial of third degree such that the first derivative 
is continuous. Apart from this detail the coupling in this model, which we will refer to 
as type B, is uniquely determined by a given value as(0) at vanishing momentum. This 
model is in general, i.e. for arbitrary 0^(0), not in agreement with the condition (|4^). 
It is not constructed to serve as a realistic description of the strong coupling. Instead, 
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Figure 3: Models for the strong coupling constant (A, B) and the behaviour according 
to one-loop renormalization (C) 



it will mainly be used in order to study the qualitative effects of Gribov's equation. 
For this purpose a model is useful in which the coupling is clearly supercritical in a 
large region of momenta. It will also be useful to vary the strength of the coupling 
and to study the effects resulting from this change. Similarly, a comparison between 
the results obtained with coupling of type A and of type B will be interesting. The 
differences will in that case also depend on the initial conditions of the solutions. 



4.3 The solutions in the Euclidean region 



We now turn to the numerical study of the pair of differential equations (p^ , ( 15 ) in 
the Euclidean region of space-like momenta. We have used two different numerical 
methods. All solutions presented below have been found using a Runge-Kutta proce- 
dure, i. e. the step-wise integration of the equation starting from a set of given initial 



conditions. We have also used the routine COLSYS [18|, a non-local collocation proce- 
dure using B-splines which is also suited for boundary- value problems with boundary 
conditions given at different points. Within numerical errors agreement has been found 
in all cases in which both methods have been applied. 

The solutions in the Euclidean region can be classified according to their behaviour 



at small and large momenta q. As discussed in section 3.2, the possible fixed points for 
g — > and q ^ oo differ from each other by the values of the function x (or equivalently 
(j)). The values of p are then fixed in the respective limits. Due to the 27r-periodicity 
of the equations in x we can restrict ourselves in the following to the case in which 
X ^ vr for — > oo. We will distinguish three classes of solutions in which the function 
X approaches in the limit q ^ the values (or 27r), —it, or vr, respectively. We will 
now discuss these three classes separately. 
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The solutions in the first class|^ start at x{q = 0) = 0. As discussed in section 
this implies p{0) = 1. The free parameter in this class of solutions is therefore 
(or, equivalently, the renormalized mass m/j, see section ^ below). Here and in the 
following we will measure q in units of GeV, and thus x' i^i units of GeV~^. For small 
x'(0) the function x approaches vr monotonically. Such a solution is shown in Fig. ^ 
and we have plotted both x ^-^d p. If the coupling constant is subcritical for all q 
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Figure 4: Solutions x and p for x'(0) = 20 (corresponding to mR = 100 MeV, see below) 
with running coupling of type A 

(for example for model B with 0^(0) < Oc) the solution is monotonic for all possible 
x'(0). This situation changes if the coupling constant is supercritical in some interval 
of momenta. For small values of x'(0) the solution is still monotonic, c. f. the solution 
in Fig. Q which is found for supercritical coupling. Due to the smallness of x'(0) these 
solutions come close to vr only at momentum scales at which the running coupling is 
already subcritical again, and we do not observe oscillations. But for larger values of 
x'(0) the function x increases more rapidly and can pass the value vr. Such a solution 
is shown in Fig. |5[ For even larger x'(0) the function x can pass vr more often, and in 
fact even arbitrarily often for x'(0) ~^ oo. However, these oscillations are very strongly 
damped as the scale on the vertical axis in Fig. ^ illustrates. In all cases the solutions 
stop oscillating as q increases. This was to be expected since the coupling becomes 
subcritical at larger momentum scales. If x'(0) is negative the function x approaches 
—TT instead of vr at large q. Due to the periodicity of the equations in x we can then 
assume that these solutions start at x(0) = ^vr (instead of 0) and approach the fixed 
point vr. Otherwise the behaviour of these solutions (oscillations, damping etc.) is not 
different from the case of positive x'(0)- 

In the second class of solutions x runs from — vr to vr as g runs from to oo. In other 
words, X goes from one minimum of the potential to a neighboring one. At the same 

^These are the (only) solutions discussed in 0]. 
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Figure 5: Solution x x'(0) = 250 (corresponding to mR = 8 MeV, see below) with 
running coupling of type A 

time p runs from + 3/3^(0) to y^l + 3/3^(oo) = 2. It is possible that the zeros of 
X and p coincide. Such a solution is presented in Fig. |^ In a sense, these solutions are 
'symmetric'. (This term would be more appropriate if the coupling was not running 
but constant.) There are also solutions in which the zeros of x P do not coincide, 
they are 'asymmetric' in this sense. 

The third class comprises such solutions in which x starts at vr and also approaches 
this value at large q. The solutions thus stay in one well of the potential. At the same 
time p runs over the same range as in the previous class of solutions. Similar to that 
case, it is possible that x = and p = coincide. Such a 'symmetric' solution is shown 
in Fig. 1^ Of course, there are also 'asymmetric' solutions in which these special values 
of X ™d p occur at different momentum scales. 

If the coupling is supercritical the solutions of the second and third class can exhibit 
oscillations of x around tt (and also around — vr for — > in the third class) similar to 
the ones described above for the first class. 

Numerically, the solution of the first class can be found by integrating the equation 
starting at g = 0. In order to integrate the solutions numerically for the second and 
third class one has to start at some intermediate %. This is no restriction since the 
matching of two solutions for < g < go and go < 9 < oo is trivial. It just reflects the 
fact that integrating a differential equation starting from a stable fixed point leads to 
exponentially large numerical errors. 

The three classes discussed above exhaust all asymptotically free solution in the 
Euclidean region, i. e. solutions which end up in one of the fixed points for q oo. 
Especially, there are no solutions in which x runs from a minimum of the potential (^) 
to any other minimum that is not a neighboring one. 
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Figure 6: Solutions x ^^'^ P foi' x(0-4) = 0, x'(0-4) = 10 and p{OA) = with running 
coupling of type A 




Figure 7: Solutions x and p for x(0-5) = tt, x'(0-5) = 4 and ^i(0.5) = with running 
coupling of type A 
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5 Mass renormalization and chiral symmetry breaking 



In this section we would like to address the question how the dynamical mass function 
M(g^) of the quark behaves in the Euclidean region for the solutions discussed in the 
preceding section. We will mainly concentrate on the first class of solutions and only 
briefly comment on the other two classes at the end of the section. 

Let us define[] the 'renormalized' mass niR as the limit of the mass function M{q^) 
as the momentum vanishes, 

ruR = lim M{q^) . (48) 
g-+0 



Since in the first class of solutions x(0) = we can expand ( [1^ ) to find that the 
renormalized mass becomes 

2 

rriR = -yj-- ■ (49) 
X (0) 

Since x'(0) just the free parameter specifying these solutions we can use tur instead 
to characterize them uniquely. Although the renormalized mass is the small-momentum 
limit of the dynamical mass function it would most probably be too simple to interpret 
it as a constituent mass of the quark. 

Further we want to define a 'perturbative' mass mp of a given solution. It is 
supposed to reflect the behaviour of the perturbative tail of the mass function of this 
solution. This could be achieved by computing its asymptotic behaviour which is 
eventually described by eq. (^). For our purposes it turns out to be sufficient and 
more convenient for the numerical study to choose a definition which involves only a 
finite momentum scale. We therefore define the perturbative mass as the value of the 
mass function at the scale A at which the coupling becomes subcritical (see also Fig. 
^) and the perturbative behaviour is expected to set in, 

mp = M(A2) = Acot^. (50) 
Obviously, the scale A in this definition depends on the model for the running coupling. 



In the models A and B discussed in section 4.2, and presumably in most other realistic 



models, the values of A are very similar. Any other choice of scale would lead to similar 
results as long as that scale is chosen larger than A. 

A general property of the solutions of Gribov's equation is the rapid decrease of the 
mass function with increasing momentum. This is illustrated in Fig. ^ for two solutions 
that are similar to the one shown in Fig. ^. Both solutions in this figure do not exhibit 
oscillations since they correspond to comparatively small values of x'(0) (see also the 
corresponding discussion in section [4.31) . 

It is now interesting to study the relation between the renormalized mass tur and 
the perturbative mass mp. First we consider the case of subcritical coupling. We 
choose model B for the running coupling with a maximal value as(0) = 0.3. Since the 
coupling is subcritical for all momenta the scale A cannot be defined via the critical 
value Oc in this model. Therefore we have to supplement the definition ( |50|) of the 

^This term was introduced in and we adopt it here. This mass is not meant to be renormahzed 
as opposed to being a bare mass in the usual sense. The mass function of course describes the mass 
renormahzation for all . This 'renormalization' is done down to g = 0, having in mind a starting 
point at large q. 
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Figure 8: The mass function M{q^) for renormalized masses mn - 
respectively, for coupling of type A 
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perturbative mass with the choice A = 1 GeV in this case. This choice is to some 
extent arbitrary. It is mainly motivated by the values of A resulting from the models 
A or B with supercritical as(0). Other choices for A in the subcritical case lead to 
similar results for the dependence of mn on mp. This dependence is shown in Fig. 
^. For subcritical coupling there is a one-to-one correspondence between renormalized 
mass mn and perturbative mass mp. The renormalized mass decreases with decreasing 
perturbative mass, and vanishes when the perturbative mass vanishes. Chiral symmetry 
is thus not broken at subcritical coupling. 

If the coupling is supercritical at small momenta some of the solutions exhibit 
oscillations. These are possible until the running coupling becomes subcritical, i. e. as 
long as < A. There are also monotonic solutions for which mp has to be sufficiently 



large (see section 4.3). Let us consider one of the latter. The corresponding function 
X has a certain value at the scale A. But there are also solutions x with smaller mp = 
2/x'(0) which oscillate and pass vr twice. It is now possible that one of these solutions 
has the same value at the scale A as the monotonic solution considered before. There 
can in fact be even more solutions with this property, among them also solutions with 
negative mp (resp. negative x'(0)). The number of these solutions will (for mp ^ 0, see 
below) remain finite. The reason for this is the following. A solution x which oscillates 
very often will be very close to vr due to the strong damping in the corresponding 
equation. As a result there is a maximal possible x(^) which the solutions can reach 
for a given number of oscillations. In summary, we find that different values of mp (resp. 
x'(0)) can lead to identical values of x(^)- But since the perturbative mass mp depends 
only on this value x(-^)! see (|50|), the correspondence between the renormalized and the 
perturbative mass is no longer one-to-one. Instead it takes the form shown in Fig. 10, 
and the renormalized mass becomes a multi-valued function of the perturbative mass. 
Here we have chosen model B for the running coupling, thus A = 1.27 GeV. Further 
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Figure 9: Renormalized mass niR vs. perturbative mass mp for subcritical coupling 
(coupling of type B with q;s(0) = 0.3) 
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Figure 10: Renormalized mass niR vs. perturbative mass mp for coupling of type B 
with as (0) = 0.94 
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we have chosen 0^(0) = 0.94 in this figure because the effect is more pronounced at 
large couphng. The dependence on the couphng strength and on the model will be 
studied below. The figure has been obtained by integrating the differential equation 
up to g = A with varying initial conditions 



Fig. 10 shows that at supercritical coupling the renormalized mass does not vanish 
in the limit of vanishing perturbative mass, and we thus find that chiral symmetry is 
broken. This figure is one of our main results and we will now study it in some more 
detail. At large perturbative mass there is only one branch of solutions, denoted as (pi 
in the figure. But if the perturbative mass is below a critical mass rric there are two 
additional solutions with different renormalized mass, denoted as the branches (/>2, </>3 in 
the figure, and we can regard this as a phase transition in the vacuum of light quarks. If 
we consider only the region of positive perturbative masses we can identify the branch 
(pi with monotonic solutions, whereas the branch (p^ corresponds to solutions in which 
X passes vr. Further, the branch (p2 can be identified with solutions in which x passes 
TT and in addition has a turning point. 



We observe that the curve in Fig. |10| is symmetric with respect to the origin. This is 
an immediate consequence of the fact that the equations (0), ( |T^ ) are invariant under 
the exchange 

(p^2m-(p. (51) 

The mass function changes sign under this transformation, and hence the symmetry. 

The phase transition at the critical value nic of the perturbative mass mp is not 
the only one. In fact there is an infinite series of similar phase transitions taking place 
in the limit mp — > 0. The curve in Fig. IC exhibits a very interesting self-similarity 



which illustrates this series of phase transitions. In Fig. 11 we show a detail of Fig. 
|Tc| around the origin. Both figures show different parts of the same curve. We denote 
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Figure 11: Renormalized mass mp vs. perturbative mass mp for coupling of type B 



with a, (0) = 0.94, detail of Fig. 10 



the critical mass of the second phase transition as m^. In the two additional solutions 
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occurring at this scale the function x passes the value tt twice. 

It is interesting to note that in each phase transition one of the two additional 
branches of solutions has the quite unusual property that the renormalized mass grows 
with decreasing perturbative mass. 

The phase transitions lead to the generation of pions as Goldstone bosons |2|, |3|. 
The Bethe-Salpeter amplitude of the pion can be shown to he tp = C{G~^ ,^5} with 
a constant C. It solves an equation for qq bound states which is derived in a similar 
approximation scheme as the equation for the Green function |^. It has been speculated 
that the observation of just one pion in nature should restrict the physical value of the 
perturbative mass of the quark to be between the first and second phase transition 
[]lS| ]. As can already be seen from Figures 10 and 11, however, two successive phase 
transitions happen to take place at mass scales which differ from each other by two 
orders of magnitude. Therefore it is not possible to deduce any considerable restriction 
on the physical perturbative mass of the quark in this way. 

It is now instructive to study the dependence of the critical mass scales on the 
value and on the model chosen for the running coupling. For this we use the two 
models A and B for the running coupling introduced in section 4.2. There we did not 
yet define how to vary the coupling strength in model A. We will use the value as(0) at 
vanishing momentum as a parameter (as in model B) and adjust the value of a in eq. 
( p7|) accordingly to obtain a running coupling of varying strength in this model as well. 
Of course, the condition ( |4^ will then no longer be fulfilled. In model B the momentum 
scale A at which the coupling becomes critical is independent of the value as(0), but 
in model A it depends by construction on the parameter a. Since this dependence is 
rather weak and in order to make the results comparable we use also for model A the 
fixed value A = 1.27 GeV which is obtained for model B. 

Fig. |l^ shows the dependence of the first critical mass rric on the coupling strength 
in models A and B. In order to obtain this figure we have computed curves similar to 
the one in Fig. 10 for different values of the coupling and determined the turning point 
corresponding to rric- Below the critical coupling there are no phase transitions and 
the critical mass acquires values different from zero only for supercritical values of the 
coupling. As expected, the critical mass then grows with increasing coupling. Realistic 
values of the coupling are roughly of the order a<j(0) ~ 0.8 (model A), corresponding 
to g{0) ~ 0.3. From Fig. |l^ we see that in this region both models give very similar 
values for the critical mass rric in the range of a few MeV. Only for perturbative quark 
masses below this small nic we expect chiral symmetry breaking. This is in perfect 
agreement with Gribov's physical picture according to which exactly the existence of 
very light quarks with masses below of few MeV leads to a dramatic change in the 
vacuum structure. 

At very large values of the coupling there is a maximum and the critical mass 
decreases again. For model B this maximum is reached at g{0) ~ 1.1, and for model 
A it is found at g{0) ~ 1.5 (outside the range shown in Fig. [T^ ). This seemingly 
strange behaviour has its origin in a property of the equation that we notices already 
in sections \i.2\ There we found that the oscillations in the solutions x around the value 
TT disappear at very large coupling g > 1 + a/2/3. Had we chosen a fixed value for as 
instead of a running Ug in Fig. [l^ we would have found in fact that the critical mass 
vanishes again at g = 1 + \/2/3. But in our models the running coupling is a continuous 
function. Even for arbitrarily large g{0) it has in some momentum range values which 
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Figure 12: Critical mass rUc of the first phase transition as a function of the couphng 
with behaviour of type A and B, respectively 



make oscillations possible, and rric does not vanish even at large coupling. Obviously, 
this is true for all possible continuous shapes of the running coupling. We therefore 
expect that the occurrence of chiral symmetry breaking due to this mechanism is largely 
independent of the details of the running coupling at small momenta. 

In Fig. H we show the second critical mass m'^ as a function of the coupling. As we 
already mentioned the values of are almost two orders of magnitude smaller than 
the corresponding values of the first critical mass rric- The behaviour of m'^ at very 
large coupling (not shown in the figure) is analogous to the behaviour observed for the 
first critical mass, and there is a similar maximum and decrease at very large ^(O). 

Finally, we would like to comment on the solutions of the second and third class 
found in section |4.3| . Here the interpretation in terms of a renormalized mass is more 
involved than in the solutions of the first class discussed above. This is because here 
we find x(0) = — tt (second class) or x(0) = vr (third class), respectively. According to 
(12) this implies that a renormalized mass defined as the limit of the mass function at 
zero momentum, see (|4^), would always vanish. Nevertheless, it should be possible to 
define a quantity similar to the renormalized mass at some intermediate but small 
momentum scale. In the second and third class of solutions the behaviour at large mo- 
menta is similar to that of the first class, especially concerning the oscillations. Given 
a suitable definition of a renormalized mass in the above sense one would therefore 
observe similar phase transitions and the breaking of chiral symmetry. A full interpre- 
tation of the renormalized mass, especially for the second and third class of solutions, 
remains to be found. Most probably this would require a better understanding of the 
difficult problem of the emergence of a constituent mass of the quark. 
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Figure 13: Critical mass m'^ of the second phase transition as a function of the couphng 
with behaviour of type A and B, respectively 




6 Analytic structure of the solutions 

So far we have discussed the properties of the solutions of Gribov's equation in the region 
of space-like momenta. Now we turn to the problem of determining the corresponding 
analytic structure of the Green function in the whole momentum plane. The locations 
and nature of the singularities in the complex plane contain crucial information about 
the solutions, especially about their properties regarding causality and confinement. It 
has been argued in |^ that a suitable solution for the retarded Green function of a 
confined quark should be free of singularities in the physical region Qmq^ > 0. The 



analysis of the first class of solutions (see section O) in [Q] has shown that we should 
not expect to find such a confining solution of the equation without pion corrections. 
We will confirm that result and extend it also to the two other classes of solutions. 

As a first step we need to define the coupling constant in the whole complex mo- 
mentum plane. This problem, however, is poorly understood even at large momenta 
since the one^oop formula for the running coupling does not possess a unique ana- 
lytic continuation. Obviously, the situation is even worse at low \q\. In order to avoid 
additional parameters and artificial singularities we will therefore make the simplest 
possible assumption and choose a fixed coupling throughout the complex q-plane. The 
value of this coupling will be chosen to be supercritical, and in the examples below 
as = 0.8. We have to hope that this choice will be a suitable approximation at least for 
small \q\. For large \q\ our choice seems less appropriate. But it will turn out that the 
singularities of the Green function occur at comparatively small values of Iq^I, probably 
indicating that an adiabatic modification of the coupling at large \q\ would not change 
the analytic structure. Obviously, the choice of a constant and supercritical coupling 
disagrees with the models for we have used to find the solutions in the Euclidean 
region. The details of the solutions will in fact change at large space-like momenta. 
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But we have checked that the classification of the solutions given in section remains 
valid. 

In order to study the analytic structure we integrate the equation along suitable 
contours in the complex plane. These should start in the Euclidean region where the 
necessary conditions on the initial values are known and the possible solutions have 
been classified (see section 14. 31). Some typical contours are shown in Fig. 14. In general 




Figure 14: Typical contours in the g-plane 

(especially for the solutions of the second and third class) the contours need not start 
in the origin of the g-plane. We should remind the reader that our definition (^) 
implies that the right-hand half of the g-plane already covers the full g^-plane, and the 



Euclidean region is given by the vertical axis in Fig. |14. 

The occurrence of poles and cuts in the analytic continuation will in general lead 
to a complicated Riemann surface. We will in the following only discuss singularities 
arising in the physical region of the q'-plane, although it is in principle possible to locate 
also further singularities on unphysical sheets of the Riemann surface. 

We start with the solutions of the first class for which (^{q = 0) = and (j) ^ in 
for large space-like q, see Fig. ^ for a typical example. The different solutions of this 
class are parametrized by the renormalized mass m^. All of these solutions have a 
pole on the real (time-like) g-axis and a cut along this axis starting at the pole. There 
are no other singularities in the physical region of the g-plan^. This result holds for 
all solutions of this class independent of their renormalized mass. In particular, this 



structure is found even for solutions which are located on different branches in Fig. 10 
The singularity and the cut on the real axis have been discussed in ||2|, |^ but there 
the absence of further singularities could not be proved. To illustrate our result we 



show in Fig. 15 the continuation of the solution of Fig. ^ i. e. in this solution again 
mR = 100 MeV. In this and in the following figures the momentum q is given in units 
of GeV. The cut is clearly visible. Below the cut the function (j) is real-valued on the 
real axis. The starting point of the cut coincides with a pole of the Green function 
at g = m* , and in our example m* = 128 MeV. The pole is most easily seen in the 
function p, the real part of which is plotted in Fig. ^ in a small region around the 
pole. At this point the function p (see eq. @) vanishes such that the Green function 
develops a pole. In general, the pole on the real axis moves to larger values of m* when 



There are further poles on the unphysical sheets of the cut which we will not discuss here. 
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Figure 15: Continuation of 9m of the solution in Fig. ^ with a cut on the real g-axis 




Figure 16: Pole of 3?ep on the real axis, continuation of the solution in Fig. ^ 
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rriR is increased. 

This analytic structure of the Green function in the first class of solutions resembles 
the analytic structure in perturbation theory. The position of the pole can be inter- 
preted as the mass m* of a propagating particle. Obviously, these solutions do not lead 
to confined quarks. 

Let us now consider the solutions of the second class in which (j) runs from —itt to itt 
in the Euclidean region. These solutions have been obtained in section [4.3| by choosing 
at some point = if on the imaginary g-axis the initial values i;^ = and p = 0. (This 
is the 'symmetric' case, for the other case see below.) We find that in all solutions of 
this kind the Green function develops a pole exactly on the circle with radius r, with 
a cut starting at the pole. The exact position of the pole depends on the value chosen 
for x'{qa)- (In this section we give this parameter of the solutions in terms of x rather 
than in terms of (p since x'{qa) is real-valued and in order to make it easier to compare 
with the preceding section. We recall that the functions are related by 4> = ix and the 
derivative is with respect to q = —iq, see eq. (^l|).) For small values of this parameter 
the pole is close to the real g-axis, and can in fact be shifted arbitrarily close to the 
real axis by choosing sufficiently small x'{qa)- For larger values of this parameter the 
pole moves along the circle towards the imaginary g-axis. In all cases the pole is found 
in the physical region above the real axis. In addition, there are two further poles and 
cuts in the physical region. These two poles are located under the same polar angle in 
the g-plane. For smaller values of x'{qa) they are located closer to the real axis, and 
the one located at smaller \q\ moves closer to the origin, whereas the other one moves 
to larger distances from the origin. At all three poles the function p vanishes. As an 
example we show in Fig. ^ the function 9m for the same solution which is shown 
in Fig. ^ in the Euclidean region, i. e. for the solution obtained from r = 0.4 GeV and 
x'{qa) = 10. The shape of the cuts in the figure is due to the choice of the contours of 




Figure 17: Qmcp in the physical region, continuation of the solution in Fig. |6| 

integration. The corresponding poles can be found again by looking at the function p. 
Their positions can be seen in Fig. ^as the points where the contour lines meet. There 
are further singularities in the unphysical region which we do not discuss here. 

There are further solutions in the second class which we called 'asymmetric' in 



section 4.3 in which the zeros of (p and p in the Euclidean region do not coincide. Also 
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these solutions have a singularity in the physical region on the circle on which p = 
is chosen on the imaginary g-axis. As in the 'symmetric' solution this pole can be 
shifted very close to the real axis but not below it. The 'asymmetric' solutions also 
have further singularities but we will not discuss them here. 

Finally, there is the third class of solutions in which </> stays close to in throughout 
the Euclidean region. Here the situation is similar to the one for the solutions of the 
second class. We find a pole and a cut on the circle on which on the imaginary g-axis 
p = 0. Again this is the case for the 'symmetric' solutions (in which p = coincides 
with (f) = in on the imaginary axis) as well as for the 'asymmetric' solutions. As an 



example we present in Fig. 18 the analytic structure of the 'symmetric' solution that we 



have shown in Fig. 0. The crosses and lines correspond to the poles and cuts. The solid 
line indicates on which sheets of the cut the other singularities are located. Denoting 
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Figure 18: Positions of the singularities in the physical region of the g-plane for the 
continuation of the solution in Fig. |^ 

again as qa the point on the imaginary axis at which (p = and p = 0, we have for this 
solution QA = i • 0.5 GeV and x'iQA) = 4. For smaller values of x'iQA) the pole moves 
towards the real axis but stays in the physical region. There are two further poles and 
cuts as indicated in the figure. These two poles are located on a straight line that goes 
through the origin. 

In summary, we have found in this section that the analytic structure of the solutions 
of the first class is similar to the analytic structure of a perturbative Green function. 
For this class we find a pole and a cut on the real g-axis and no other singularities in 
the physical region. In the other two classes there are poles and cuts on the physical 
sheet. Those solutions will in general not permit an interpretation in agreement with 
causality and unitarity and should therefore be regarded as unphysical solutions of 
Gribov's equation. We did not find any hint for the existence of exceptional solutions 
in which those singularities move to the unphysical sheet for special initial conditions. 
We can thus conclude that, as anticipated in |^], there are no solutions of the equation 
without pion corrections which would lead to confined quarks. 

7 Summary and outlook 

The subject of this paper has been Gribov's equation for the Green function of light 
quarks. We have have outlined how this equation is derived from the Dyson-Schwinger 
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equation, and we have described how the approximations made in the derivation are 
motivated by the physical picture of supercritical color charges in QCD. In Gribov's 
scenario the phenomenon of supercritical charges can occur in QCD due to the existence 
of very light quarks and is expected to cause chiral symmetry breaking and confinement. 
In the present paper we have concentrated on the equation proposed in ||2|, ^ which 
does not yet contain pion corrections. The equation describes the Green function in 
Feynman gauge. It collects the most singular contributions to the Dyson-Schwinger 
equation coming from the region of small momenta. At the same time it reproduces 
asymptotic freedom at large momenta. The derivation of the equation involves the 
assumption that the running coupling remains finite at small momentum scales. This is 
in agreement with recent phenomenological results obtained in the dispersive approach 
to power corrections in QCD. These analyses find a typical value of the coupling in the 
infrared which is above the critical value Oc = 0.43 arising in Gribov's scenario. 

Gribov's differential equation for the light quark's Green function is a nonlinear 
equation and so far only asymptotic expansions had been used to study it. We have 
described the corresponding results, especially the fixed point structure of the equa- 
tion resulting from this analysis and the emergence of the critical value of the strong 
coupling. We have then performed a complete numerical study of the equation. After 
defining suitable models for the behaviour of the strong coupling constant in the in- 
frared region we have classified the possible solutions of the equation according to their 
behaviour in the Euclidean region of space-like momenta. At supercritical coupling the 
solutions change their behaviour and can oscillate while approaching the fixed points. 
The dynamical mass function of the quark has been computed in order to show that 
this change leads to the breaking of chiral symmetry at supercritical coupling. Chiral 
symmetry breaking takes place independently of the details of the running coupling 
in the infrared as long as it is supercritical in some interval of the momentum. The 
breaking of chiral symmetry is connected with the occurrence of a series of phase tran- 
sitions in the vacuum of light quarks. If the perturbative mass mp of the quark, i.e. 
the mass defined at a large momentum scale, is below a critical mass then there are 
several solutions leading to different quark masses at low momentum scales, and in 
general rriR does not vanish even in the limit of vanishing perturbative mass mp. We 
have determined the critical mass as a function of the coupling. It turns out to depend 
considerably on the mean value of the coupling in the infrared, but the dependence is 
weaker if the coupling is only slightly above the critical value. In agreement with the 
physical picture of supercritical color charges the phase transitions occur only for very 
small values of the perturbative quark mass. 

In ||3|, ^ it has been advocated that a confining Green function in Gribov's picture 
should be free of singularities in the physical region including time-like momenta. We 
have therefore studied the solutions of Gribov's equation in the whole complex mo- 
mentum plane. One class of solutions has an analytic structure similar to that of a 
perturbative Green function, whereas the other two classes of solutions have singulari- 
ties on the physical sheet. As anticipated in Q there are no solutions of the equation 
studied in the present paper exhibiting the structure required for confinement. 

The phase transitions connected with the breaking of chiral symmetry lead to the 
generation of pions as Goldstone bosons. Their special properties in this approach Q 
indicate that they should be taken into account as elementary objects, and that the 
equation for the Green function of the light quark should be modified accordingly. This 
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modification is expected to change the behaviour of the solutions in the Euchdean 
region only very little. In particular, the breaking of chiral symmetry will take place 
in a very similar way. But the pion corrections to the equation are expected to have 
significant effects on the analytic structure of the solutions. Due to that it appears pos- 
sible to find a confining Green function. The next step should therefore be a numerical 
study of the corresponding differential equation. Although the modified equation with 
pion corrections is more complicated it can be analysed using the same methods that 
have been used in the present paper. It would at a later stage also be interesting to in- 
vestigate the effects of an additional scalar (cr-)meson which has been widely discussed 
in the literature in the context of meson spectra and of sigma models. 

We hope that the results of our analysis will be useful also for the more conventional 
approach to the Dyson-Schwinger equation for the quark. Some general properties of 
the Green function resulting from the Dyson-Schwinger equation appear to be rather 
universal and largely independent of the particular approximation scheme that is used. 
The analysis of the Dyson-Schwinger equation in quenched supercritical QED ||2^ , for 
example, leads to a picture of phase transitions and chiral symmetry breaking that 
is very similar to the one resulting from Gribov's equation. It would certainly be 
useful to identify such universal features of the Green function and to achieve a better 
understanding of their origin. In this respect it is very important to consider different 
approximation schemes among which Gribov's approach is certainly exceptional since 
it is motivated by the very interesting physical picture of supercritical charges. 
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